Protein subcellular relocalization and function of duplicated flagellar calcium binding protein genes in honey bee trypanosomatid parasite.

The honey bee trypanosomatid parasite, Lotmaria passim, contains two genes that encode the flagellar calcium binding protein (FCaBP) through tandem duplication in its genome. FCaBPs localize in the flagellum and entire body membrane of L. passim through specific N-terminal sorting sequences. This finding suggests that this is an example of protein subcellular relocalization resulting from gene duplication, altering the intracellular localization of FCaBP. However, this phenomenon may not have occurred in Leishmania, as one or both of the duplicated genes have become pseudogenes. Multiple copies of the FCaBP gene are present in several Trypanosoma species and Leptomonas pyrrhocoris, indicating rapid evolution of this gene in trypanosomatid parasites. The N-terminal flagellar sorting sequence of L. passim FCaBP1 is in close proximity to the BBSome complex, while that of Trypanosoma brucei FCaBP does not direct GFP to the flagellum in L. passim. Deletion of the two FCaBP genes in L. passim affected growth and impaired flagellar morphogenesis and motility, but it did not impact host infection. Therefore, FCaBP represents a duplicated gene with a rapid evolutionary history that is essential for flagellar structure and function in a trypanosomatid parasite.

Aspartyl protease in the secretome of honey bee trypanosomatid parasite contributes to infection of bees.

BACKGROUND: The exoproteome, which consists of both secreted proteins and those originating from cell surfaces and lysed cells, is a critical component of trypanosomatid parasites, facilitating interactions with host cells and gut microbiota. However, its specific roles in the insect hosts of these parasites remain poorly understood. METHODS: We conducted a comprehensive characterization of the exoproteome in Lotmaria passim, a trypanosomatid parasite infecting honey bees, under culture conditions. We further investigated the functions of two conventionally secreted proteins, aspartyl protease (LpAsp) and chitinase (LpCht), as representative models to elucidate the role of the secretome in L. passim infection of honey bees. RESULTS: Approximately 48% of L. passim exoproteome proteins were found to share homologs with those found in seven Leishmania spp., suggesting the existence of a core exoproteome with conserved functions in the Leishmaniinae lineage. Bioinformatics analyses suggested that the L. passim exoproteome may play a pivotal role in interactions with both the host and its microbiota. Notably, the deletion of genes encoding two secretome proteins revealed the important role of LpAsp, but not LpCht, in L. passim development under culture conditions and its efficiency in infecting the honey bee gut. CONCLUSIONS: Our results highlight the exoproteome as a valuable resource for unraveling the mechanisms employed by trypanosomatid parasites to infect insect hosts by interacting with the gut environment.

Evolutionary history of metazoan TMEM16 family.

Most of Transmembrane protein 16 (TMEM16) family members function as either a Ca2+-activated Cl- channel (CaCC) or phospholipid scramblase (CaPLSase) and play diverse physiological roles. It is well conserved in eukaryotes; however, the origin and evolution of different subfamilies in Metazoa are not yet understood. To uncover the evolutionary history of the TMEM16 family, we analyzed 398 proteins from 74 invertebrate species using evolutionary genomics. We found that the TMEM16C-F and J subfamilies are vertebrate-specific, but the TMEM16A/B, G, H, and K subfamilies are ancient and present in many, but not all metazoan species. The most ancient subfamilies in Metazoa, TMEM16L and M, are only maintained in limited species. TMEM16N and O are Cnidaria- and Ecdysozoa-specific subfamilies, respectively, and Ctenophora, Xenacoelomorpha, and Rotifera contain species-specific proteins. We also identified TMEM16 genes that are closely linked together in the genome, suggesting that they have been generated via recent gene duplication. The anoctamin domain structures of invertebrate-specific TMEM16 proteins predicted by AlphaFold2 contain conserved Ca2+-binding motifs and permeation pathways with either narrow or wide inner gates. The inner gate distance of TMEM16 protein may have frequently switched during metazoan evolution, and thus determined the function of the protein as either CaCC or CaPLSase. These results demonstrate that TMEM16 family has evolved by gene gain and loss in metazoans, and the genes have been generally under purifying selection to maintain protein structures and physiological functions.

DWV 3C Protease Uncovers the Diverse Catalytic Triad in Insect RNA Viruses.

Deformed wing virus (DWV) is the most prevalent Iflavirus that is infecting honey bees worldwide. However, the mechanisms of its infection and replication in host cells are poorly understood. In this study, we analyzed the structure and function of DWV 3C protease (3Cpro), which is necessary for the cleavage of the polyprotein to synthesize mature viral proteins. Thus, it is one of the nonstructural viral proteins essential for the replication. We found that the 3Cpros of DWV and picornaviruses share common enzymatic properties, including sensitivity to the same inhibitors, such as rupintrivir. The predicted structure of DWV 3Cpro by AlphaFold2, the predicted rupintrivir binding domain, and the protease activities of mutant proteins revealed that it has a Cys-His-Asn catalytic triad. Moreover, 3Cpros of other Iflaviruses and Dicistrovirus appear to contain Asn, Ser, Asp, or Glu as the third residue of the catalytic triad, suggesting diversity in insect RNA viruses. Both precursor 3Cpro with RNA-dependent RNA polymerase and mature 3Cpro are present in DWV-infected cells, suggesting that they may have different enzymatic properties and functions. DWV 3Cpro is the first 3Cpro characterized among insect RNA viruses, and our study uncovered both the common and unique characteristics among 3Cpros of Picornavirales. Furthermore, it would be possible to use the specific inhibitors of DWV 3Cpro to control DWV infection in honey bees in future. IMPORTANCE The number of managed honey bee (Apis mellifera) colonies has considerably declined in many developed countries in the recent years. Deformed wing virus (DWV) vectored by the mites is the major threat to honey bee colonies and health. To give insight into the mechanism of DWV replication in the host cells, we studied the structure-function relationship of 3C protease (3Cpro), which is necessary to cleave a viral polyprotein at the specific sites to produce the mature proteins. We found that the overall structure, some inhibitors, and processing of 3Cpro are shared between Picornavirales; however, there is diversity in the catalytic triad. DWV 3Cpro is the first viral protease characterized among insect RNA viruses and reveals the evolutionary history of 3Cpro among Picornavirales. Furthermore, DWV 3Cpro inhibitors identified in our study could also be applied to control DWV in honey bees in future.

DWV Infection in vitro Using Honey Bee Pupal Tissue.

The deformed wing virus (DWV) has been best characterized among honey bee viruses; however, very little is known regarding the mechanisms of viral infection and replication due to the lack of immortalized honey bee cell lines. To solve this problem, we established an in vitro system using honey bee pupal tissue to reconstruct DWV binding and entry into the host cell, followed by translation of the RNA genome and polyprotein processing using RNA-dependent RNA polymerase (RdRP) as a marker. Using this system, the P-domain of the virion subunit VP1 was found to be essential for DWV infection, but not for binding and entry into the cell. DWV efficiently infected the head tissue derived from early but not late pupa, suggesting that undifferentiated cells are targeted for viral infection. Furthermore, we found that inhibitors of mammalian picornavirus 3C-protease, rupintrivir and quercetin suppressed RdRP synthesis, indicating that this in vitro system is also useful for screening a compound to control viral infection. Our in vitro system may help to understand the mechanism of DWV infection in host cells.

Honey Bee Suppresses the Parasitic Mite Vitellogenin by Antimicrobial Peptide.

The negative effects of honey bee parasitic mites and deformed wing virus (DWV) on honey bee and colony health have been well characterized. However, the relationship between DWV and mites, particularly viral replication inside the mites, remains unclear. Furthermore, the physiological outcomes of honey bee immune responses stimulated by DWV and the mite to the host (honey bee) and perhaps the pathogen/parasite (DWV/mite) are not yet understood. To answer these questions, we studied the tripartite interactions between the honey bee, Tropilaelaps mercedesae, and DWV as the model. T. mercedesae functioned as a vector for DWV without supporting active viral replication. Thus, DWV negligibly affected mite fitness. Mite infestation induced mRNA expression of antimicrobial peptides (AMPs), Defensin-1 and Hymenoptaecin, which correlated with DWV copy number in honey bee pupae and mite feeding, respectively. Feeding T. mercedesae with fruit fly S2 cells heterologously expressing honey bee Hymenoptaecin significantly downregulated mite Vitellogenin expression, indicating that the honey bee AMP manipulates mite reproduction upon feeding on bee. Our results provide insights into the mechanism of DWV transmission by the honey bee parasitic mite to the host, and the novel role of AMP in defending against mite infestation.

Trypanosomatid parasite dynamically changes the transcriptome during infection and modifies honey bee physiology.

It is still not understood how honey bee parasite changes the gene expression to adapt to the host environment and how the host simultaneously responds to the parasite infection by modifying its own gene expression. To address this question, we studied a trypanosomatid, Lotmaria passim, which can be cultured in medium and inhabit the honey bee hindgut. We found that L. passim decreases mRNAs associated with protein translation, glycolysis, detoxification of radical oxygen species, and kinetoplast respiratory chain to adapt to the anaerobic and nutritionally poor honey bee hindgut during the infection. After the long term infection, the host appears to be in poor nutritional status, indicated by the increase and decrease of take-out and vitellogenin mRNAs, respectively. Simultaneous gene expression profiling of L. passim and honey bee during infection by dual RNA-seq provided insight into how both parasite and host modify their gene expressions.

Gene Disruption of Honey Bee Trypanosomatid Parasite, Lotmaria passim, by CRISPR/Cas9 System.

Two trypanosomatid species, Lotmaria passim and Crithidia mellificae, have been shown to parasitize honey bees to date. L. passim appears to be more prevalent than C. mellificae and specifically infects the honey bee hindgut. Although the genomic DNA has been sequenced, the effects of infection on honey bee health and colony are poorly understood. To identify the genes that are important for infecting honey bees and to understand their functions, we applied the CRISPR/Cas9 system to establish a method to manipulate L. passim genes. By electroporation of plasmid DNA and subsequent selection by drug, we first established an L. passim clone expressing tdTomato or Cas9. We also successfully disrupted the endogenous miltefosine transporter and tyrosine aminotransferase genes by replacement with drug (hygromycin) resistant gene using the CRISPR/Cas9-induced homology-directed repair pathway. The L. passim clone expressing fluorescent marker, as well as the simple method for editing specific genes, could become useful approaches to understand the underlying mechanisms of honey bee-trypanosomatid parasite interactions.

Honey Bee Parasitic Mite Contains the Sensilla-Rich Sensory Organ on the Foreleg Tarsus Expressing Ionotropic Receptors With Conserved Functions.

Honey bee parasitic mites (Tropilaelaps mercedesae and Varroa destructor) detect temperature, humidity, and odor but the underlying sensory mechanisms are poorly understood. To uncover how T. mercedesae responds to environmental stimuli inside a hive, we first identified the sensilla-rich sensory organ on the foreleg tarsus. The organ appeared to correspond to Haller's organ in ticks and contained four types of sensilla, which may respond to different stimuli based on their morphology. We searched for differentially expressed genes between the forelegs and hindlegs to identify mRNAs potentially associated with the sensory organ. The forelegs were enriched with mRNAs encoding sensory proteins such as ionotropic receptors (IRs) and gustatory receptors, as well as proteins involved in ciliary transport. We also found that T. mercedesae IR25a and IR93a were capable of rescuing temperature and humidity preference defects in Drosophila melanogaster IR25a and IR93a mutants. These results demonstrate that the structures and physiological functions of ancient IRs have been conserved during arthropod evolution. Our study provides insight into the sensory mechanisms of honey bee parasitic mites, as well as potential targets for methods to control the most serious honey bee pest.

Genomes of trombidid mites reveal novel predicted allergens and laterally transferred genes associated with secondary metabolism.

Background: Trombidid mites have a unique life cycle in which only the larval stage is ectoparasitic. In the superfamily Trombiculoidea ("chiggers"), the larvae feed preferentially on vertebrates, including humans. Species in the genus Leptotrombidium are vectors of a potentially fatal bacterial infection, scrub typhus, that affects 1 million people annually. Moreover, chiggers can cause pruritic dermatitis (trombiculiasis) in humans and domesticated animals. In the Trombidioidea (velvet mites), the larvae feed on other arthropods and are potential biological control agents for agricultural pests. Here, we present the first trombidid mites genomes, obtained both for a chigger, Leptotrombidium deliense, and for a velvet mite, Dinothrombium tinctorium. Results: Sequencing was performed using Illumina technology. A 180 Mb draft assembly for D. tinctorium was generated from two paired-end and one mate-pair library using a single adult specimen. For L. deliense, a lower-coverage draft assembly (117 Mb) was obtained using pooled, engorged larvae with a single paired-end library. Remarkably, both genomes exhibited evidence of ancient lateral gene transfer from soil-derived bacteria or fungi. The transferred genes confer functions that are rare in animals, including terpene and carotenoid synthesis. Thirty-seven allergenic protein families were predicted in the L. deliense genome, of which nine were unique. Preliminary proteomic analyses identified several of these putative allergens in larvae. Conclusions: Trombidid mite genomes appear to be more dynamic than those of other acariform mites. A priority for future research is to determine the biological function of terpene synthesis in this taxon and its potential for exploitation in disease control.

HsTRPA of the Red Imported Fire Ant, Solenopsis invicta, Functions as a Nocisensor and Uncovers the Evolutionary Plasticity of HsTRPA Channels.

Solenopsis invicta, the red imported fire ant, represents one of the most devastating invasive species. To understand their sensory physiology, we identified and characterized their Hymenoptera-specific (Hs) TRPA channel, SiHsTRPA. Consistent with the sensory functions of SiHsTRPA, it is activated by heat, an electrophile, and an insect repellent. Nevertheless, SiHsTRPA does not respond to most of the honey bee ortholog (AmHsTRPA)-activating compounds. The jewel wasp ortholog (NvHsTRPA) is activated by these compounds even though it outgroups both AmHsTRPA and SiHsTRPA. Characterization of AmHsTRPA/SiHsTRPA chimeric channels revealed that the amino acids in the N terminus, as well as ankyrin repeat 2 (AR2) of AmHsTRPA, are essential for the response to camphor. Furthermore, amino acids in ARs 3 and 5-7 were specifically required for the response to diallyl disulfide. Thus, amino acid substitutions in the corresponding domains of SiHsTRPA during evolution would be responsible for the loss of chemical sensitivity. SiHsTRPA-activating compounds repel red imported fire ants, suggesting that SiHsTRPA functions as a sensor for noxious compounds. SiHsTRPA represents an example of the species-specific modulation of orthologous TRPA channel properties by amino acid substitutions in multiple domains, and SiHsTRPA-activating compounds could be used to develop a method for controlling red imported fire ants.

Characterization of the Copy Number and Variants of Deformed Wing Virus (DWV) in the Pairs of Honey Bee Pupa and Infesting Varroa destructor or Tropilaelaps mercedesae.

Recent honey bee colony losses, particularly during the winter, have been shown to be associated with the presence of both ectoparasitic mites and Deformed Wing Virus (DWV). Whilst the role of Varroa destructor mites as a viral vector is well established, the role of Tropilaelaps mercedesae mites in viral transmission has not been fully investigated. In this study, we tested the effects that V. destructor and T. mercedesae infestation have on fluctuation of the DWV copy number and alteration of the virus variants in honey bees by characterizing individual pupae and their infesting mites. We observed that both mite species were associated with increased viral copy number in honey bee pupae. We found a positive correlation between DWV copy number in pupae and copy number in infesting mites, and the same DWV type A variant was present in either low or high copy number in both honey bee pupae and infesting V. destructor. These data also suggest that variant diversity is similar between honey bee pupae and the mites that infest them. These results support a previously proposed hypothesis that DWV suppresses the honey bee immune system when virus copy number reaches a specific threshold, promoting greater replication.

Draft genome of the honey bee ectoparasitic mite, Tropilaelaps mercedesae, is shaped by the parasitic life history.

The number of managed honey bee colonies has considerably decreased in many developed countries in recent years and ectoparasitic mites are considered as major threats to honey bee colonies and health. However, their general biology remains poorly understood. We sequenced the genome of Tropilaelaps mercedesae, the prevalent ectoparasitic mite infesting honey bees in Asia, and predicted 15 190 protein-coding genes that were well supported by the mite transcriptomes and proteomic data. Although amino acid substitutions have been accelerated within the conserved core genes of two mites, T. mercedesae and Metaseiulus occidentalis, T. mercedesae has undergone the least gene family expansion and contraction between the seven arthropods we tested. The number of sensory system genes has been dramatically reduced, but T. mercedesae contains all gene sets required to detoxify xenobiotics. T. mercedesae is closely associated with a symbiotic bacterium (Rickettsiella grylli-like) and Deformed Wing Virus, the most prevalent honey bee virus. T. mercedesae has a very specialized life history and habitat as the ectoparasitic mite strictly depends on the honey bee inside a stable colony. Thus, comparison of the genome and transcriptome sequences with those of a tick and free-living mites has revealed the specific features of the genome shaped by interaction with the honey bee and colony environment. Genome and transcriptome sequences of T. mercedesae, as well as Varroa destructor (another globally prevalent ectoparasitic mite of honey bee), not only provide insights into the mite biology, but may also help to develop measures to control the most serious pests of the honey bee.

TRPA1 Channels in Drosophila and Honey Bee Ectoparasitic Mites Share Heat Sensitivity and Temperature-Related Physiological Functions.

The transient receptor potential cation channel, subfamily A, member 1 (TRPA1) is conserved between many arthropods, and in some has been shown to function as a chemosensor for noxious compounds. Activation of arthropod TRPA1 channels by temperature fluctuations has been tested in only a few insect species, and all of them were shown to be activated by heat. The recent identification of chemosensitive TRPA1 channels from two honey bee ectoparasitic mite species (VdTRPA1 and TmTRPA1) have provided an opportunity to study the temperature-dependent activation and the temperature-associated physiological functions of TRPA1 channels in non-insect arthropods. We found that both mite TRPA1 channels are heat sensitive and capable of rescuing the temperature-related behavioral defects of a Drosophila melanogaster trpA1 mutant. These results suggest that heat-sensitivity of TRPA1 could be conserved between many arthropods despite its amino acid sequence diversity. Nevertheless, the ankyrin repeats (ARs) 6 and 7 are well-conserved between six heat-sensitive arthropod TRPA1 channels and have critical roles for the heat activation of VdTRPA1.

Isoform-specific modulation of the chemical sensitivity of conserved TRPA1 channel in the major honeybee ectoparasitic mite, Tropilaelaps mercedesae.

We identified and characterized the TRPA1 channel of Tropilaelaps mercedesae (TmTRPA1), one of two major species of honeybee ectoparasitic mite. Three TmTRPA1 isoforms with unique N-terminal sequences were activated by heat, and the isoform highly expressed in the mite's front legs, TmTRPA1b, was also activated by 27 plant-derived compounds including electrophiles. This suggests that the heat- and electrophile-dependent gating mechanisms as nocisensitive TRPA1 channel are well conserved between arthropod species. Intriguingly, one TmTRPA1 isoform, TmTRPA1a, was activated by only six compounds compared with two other isoforms, demonstrating that the N-terminal sequences are critical determinants for the chemical sensitivity. This is the first example of isoform-specific modulation of chemical sensitivity of TRPA1 channel in one species. α-terpineol showed repellent activity towards T. mercedesae in a laboratory assay and repressed T. mercedesae entry for reproduction into the brood cells with fifth instar larvae in hives. Thus, α-terpineol could be used as the potential compound to control two major honeybee ectoparasitic mites, T. mercedesae and Varroa destructor, in the apiculture industry.

Plant-Derived Tick Repellents Activate the Honey Bee Ectoparasitic Mite TRPA1.

We have identified and characterized the TRPA1 channel of Varroa destructor (VdTRPA1), a major ectoparasitic mite of honey bee. One of the two VdTRPA1 isoforms, VdTRPA1L, was activated by a variety of plant-derived compounds, including electrophilic compounds, suggesting that chemical activation profiles are mostly shared between arthropod TRPA1 channels. Nevertheless, carvacrol and α-terpineol activated VdTRPA1L but not a honey bee noxious-stimuli-sensitive TRPA, AmHsTRPA, and Drosophila melanogaster TRPA1. Activation of VdTRPA1L in D. melanogaster taste neurons by the above compounds was sufficient to modify the gustatory behaviors. Carvacrol and α-terpineol repelled V. destructor in a laboratory assay, and α-terpineol repressed V. destructor entry for reproduction into the brood cells in hives. Understanding the functions of parasite TRP channels not only gives clues about the evolving molecular and cellular mechanisms of parasitism but also helps in the development of control methods.

Differential diagnosis of the honey bee trypanosomatids Crithidia mellificae and Lotmaria passim.

Trypanosomatids infecting honey bees have been poorly studied with molecular methods until recently. After the description of Crithidia mellificae (Langridge and McGhee, 1967) it took about forty years until molecular data for honey bee trypanosomatids became available and were used to identify and describe a new trypanosomatid species from honey bees, Lotmaria passim (Evans and Schwarz, 2014). However, an easy method to distinguish them without sequencing is not yet available. Research on the related bumble bee parasites Crithidia bombi and Crithidia expoeki revealed a fragment length polymorphism in the internal transcribed spacer 1 (ITS1), which enabled species discrimination. In search of fragment length polymorphisms for differential diagnostics in honey bee trypanosomatids, we studied honey bee trypanosomatid cell cultures of C. mellificae and L. passim. This research resulted in the identification of fragment length polymorphisms in ITS1 and ITS1-2 markers, which enabled us to develop a diagnostic method to differentiate both honey bee trypanosomatid species without the need for sequencing. However, the amplification success of the ITS1 marker depends probably on the trypanosomatid infection level. Further investigation confirmed that L. passim is the dominant species in Belgium, Japan and Switzerland. We found C. mellificae only rarely in Belgian honey bee samples, but not in honey bee samples from other countries. C. mellificae was also detected in mason bees (Osmia bicornis and Osmia cornuta) besides in honey bees. Further, the characterization and comparison of additional markers from L. passim strain SF (published as C. mellificae strain SF) and a Belgian honey bee sample revealed very low divergence in the 18S rRNA, ITS1-2, 28S rRNA and cytochrome b sequences. Nevertheless, a variable stretch was observed in the gp63 virulence factor.

Evolutionary dynamics of metazoan TRP channels.

Transient receptor potential (TRP) channels are unusual among cation channels because of their diverse cation selectivities and activation mechanisms. TRP channels thus play major roles in various sensory perceptions by functioning as multimodal signal integrators. Some TRP subfamily members are also implicated in acute and chronic pain and inflammation. So far, most TRP channel studies have been targeted to human and model organisms within a limited evolutionary context. Classification of TRP channels in various animal genomes has revealed extensive gene gain and loss events across animal species. Furthermore, the chemical activation profiles of some orthologous TRP channels were different between species such as human and mouse. Amino acid substitutions must underlie such differences, and the crucial amino acid residues have been identified in some cases. These changes represent the evolution of TRP channels at the amino acid sequence level. There is also evidence that TRP channels have obtained species-diversity through alternative splicing and possibly cis-regulatory element mutations. All of the above demonstrate the dynamic and plastic evolutionary history of metazoan TRP channels at multiple levels, possibly in conjunction with the specific habitats and life histories of individual species.

Evolution of TRP channels inferred by their classification in diverse animal species.

The functions of TRP channels have primarily been characterized in model organisms within a limited evolutionary context. We thus characterize the TRP channels in choanoflagellate, sponge, Cnidaria, Lophotrochozoa, and arthropods to understand how they emerged during early evolution of animals and have changed during diversification of various species. As previously reported, five metazoan TRP subfamily members (TRPA, TRPC, TRPM, TRPML, and TRPV) were identified in choanoflagellates, demonstrating that they evolved before the emergence of multicellular animals. TRPN was identified in Hydra magnipapillata, and therefore emerged in the last common ancestor of Cnidaria-Bilateria. A novel subfamily member (TRPVL) was identified in Cnidaria and Capitella teleta, indicating that it was present in the last common ancestor of Cnidaria-Bilateria but has since been lost in most bilaterians. The characterization of arthropod TRP channels revealed that Daphnia pulex and insects have specifically expanded the TRPA subfamily, which diverged from the ancient TRPA1 channel gene. The diversity of TRPA channels except TRPA1 was detectable even within a single insect family, namely the ant lineage. The present study demonstrates the evolutionary history of TRP channel genes, which may have diverged in conjunction with the specific habitats and life histories of individual species.

Molecular detection of protozoan parasites infecting Apis mellifera colonies in Japan.

The role of protozoan parasites in honey bee health and distribution in the world is not well understood. Therefore, we carried out a molecular survey for the presence of Crithidia mellificae and Apicystis bombi in the colonies of both non-native Apis mellifera and native Apis cerana japonica in Japan. We found that A. mellifera, but not A. c. japonica, colonies are parasitized with C. mellificae and A. bombi. Their absence in A. c. japonica colonies indicates that A. mellifera is their native host. Nevertheless, the prevalence in A. mellifera colonies is low compared with other pathogens such as viruses and Nosema microsporidia. Japanese C. mellificae isolates share well-conserved nuclear-encoded gene sequences with Swiss and US isolates. We have found two Japanese haplotypes (A and B) with two nucleotide differences in the kinetoplast-encoded cytochrome b sequence. The haplotype A is identical to Swiss isolate. These results demonstrate that C. mellificae and A. bombi distribute in Asia, Oceania, Europe, and South and North Americas.